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ABSTRACT 

Massive black holes in galactic nuclei vary their mass Mbh and spin vector Jbh due to accretion. In 
this study we relax, for the first time, the assumption that accretion can be either chaotic, i.e. when 
the accretion episodes are randomly and isotropically oriented, or coherent, i.e. when they occur all in 
a preferred plane. Instead, we consider different degrees of anisotropy in the fueling, never confining 
to accretion events on a fixed direction. We follow the black hole growth evolving contemporarily 
mass, spin modulus a and spin direction. We discover the occurrence of two regimes. An early 
phase (Mbh ^ 10^ M©) in which rapid alignment of the black hole spin direction to the disk angular 
momentum in each single episode leads to erratic changes in the black hole spin orientation and at 
the same time to large spins (a ^ 0.8). A second phase starts when the black hole mass increases 
above > 10^ M© and the accretion disks carry less mass and angular momentum relatively to the hole. 
In the absence of a preferential direction the black holes tend to spin-down in this phase. However, 
when a modest degree of anisotropy in the fueling process (still far from being coherent) is present, 
the black hole spin can increase up to a ^ 1 for very massive black holes (Mbh ^ 10^ Mq), and its 
direction is stable over the many accretion cycles. We discuss the implications that our results have 
in the realm of the observations of black hole spin and jet orientations. 

Subject headings: Black hole physics — Accretion, accretion disks — galaxies: active — galaxies: jets 



1. INTRODUCTION 

Massive black holes (BHs) that inhabit in the nuclei of 
massive galaxies are described by only two parameters: 
their mass Mbh and spin Jbh- 

Tens of BH masses in quiescent, nearby galaxies have 
been directly estimated to date studying the dynam- 
ics of stars and gas and these masses all fall within 
> IO^Mq - 101°Mq (e.g. Giiltekin et al. 2009, Mc- 
Connell et al. 2011 and references therein). Similarly, 
tens of thousands BH masses in active nuclei have also 
been estimated using empirical relations, and, while for 
quasars their median mass is in excess of 10^ M© (e.g., 
Vestergaard et al. 2008), faint AGN are powered by BHs 
with masses as small as ^ 10^ M© (Peterson et al. 2005). 

Measurements of the dimensionless spin parameter a 
are more controversial. The parameter a has been re- 
cently measured only for few AGNs (e.g. Brenneman & 
Reynolds 2006; Schmoh et al. 2009; de la Calle Perez et 
al. 2010; Patrick et al. 2011a,b; Gallo et al. 2011; Bren- 
neman et al. 2011) through X-ray spectroscopy. How- 
ever, several free parameters enter in the estimate, creat- 
ing a severe degeneracy problem. For example, different 
groups obtained different esitmates of a in MCG-6-30-15: 
a > 0.98 (Brenneman & Reynolds 2006), a = 0.86 ±0.01 
(de la Calle Perez et al. 2010) and a = 0.491°;^° (Patrick 
et al. 2011a). Even more problematic is the case of NGC 
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3783, for which different groups found a > 0.88 (Bren- 
neman et al. 2011) and a < 0.32 (Patrick et al. 2011b). 
Measuring and understanding spins is crucial to assess 
the cosmic evolution of massive BHs. Firstly, spins af- 
fect the accretion-luminosity conversion efficiency; highly 
spinning BHs can convert up to ^ 40% of the accreted 
matter into radiation, making them more luminous, al- 
beit making their growth slowlier. Secondly, the spin 
paradigm assumes that radio jets observed in AGNs are 
launched by highly spinning BHs (Blandford & Znajek 
1977). Lastly, spins dramatically afSfect the gravitational 
recoil suffered by the remnant BH after a binary merger. 
It has in fact been shown that highly spinning BHs can 
experience kicks up to 5000 km s~^ depending on their 
progenitor spin magnitude and orientation (Campanelli 
et al. 2007; Baker et al. 2008; Herrmann et al. 2007; 
Schnittman & Buonanno 2007; Lousto & Zlochower 2011; 
Lousto et al. 2012). These super-kicks are sufficient 
to eject the remnant from the deepest potential well of 
the most massive galaxies, with potentially important 
implications for the occupation fraction of massive BHs 
in galaxies (Schnittman 2007; Volonteri 2007; Volonteri, 
Haardt & Giihekin 2008; Volonteri, Giihekin & Dotti 
2010). 

BH masses and spins are thought to build up through 
gas accretion and BH mergers, and their history of 
growth erases the initial values of Mbh and a. The main 
driver of the BH spin evolution is gas accretion (Berti & 
Volonteri 2008; Fanidakis et al. 2011; Barausse 2012)0. 
Two different accretion scenarios have been proposed to 
date: (i) coherent accretion, in which the BHs accrete gas 
with a well defined, almost constant, angular momentum 
direction, and (ii) chaotic accretion, in which parcels of 
gas accrete on the BHs in randomly oriented planes (e.g. 

^ With the possible exception of the most massive BHs in mas- 
sive, gas-poor, low-redshift ellipticals (e.g. Fanidakis et al. 2011). 
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King et al. 2005; King & Pringlc 2006; 2007). The two 
models result in different BH evolutions and different 
expected distributions of the spin magnitudes. Coher- 
ent accretion keeps on adding angular momentum to the 
BHs in the same direction and results in very high spins 
(0.8 < a < 1) aligned with the angular momentum of the 
accreting material (e.g. Dotti ct al. 2010 and references 
therein). The chaotic case is more subtle: gas accretion 
on a rotating BH on a retrograde orbit has a larger last 
stable orbit than gas with angular momentum aligned 
with the BH spin. As a consequence, retrograde accret- 
ing gas transfers to the BH more (and negative) angular 
momentum per unit of mass than the progradc. If ret- 
rograde and progradc accretion arc equally probable (as 
implicitly assumed in the chaotic scenario), the BH spin 
a is biased toward low values (0 < a < 0.2, King, Pringle 
& Hofmann 2008; Berti & Volonteri 2008). 

These two models cover the extreme cases in which the 
gas either flows from a stable, fixed direction or from fully 
random directions. In this investigation we relax these 
extreme and unrealistic assumptions, exploring the evo- 
lution of the BH mass and spin vector varying the degree 
of anisotropy in the fueling process to mimic anisotropics 
present in the gas in the nuclear regions of active galax- 
ies. We demonstrate how the evolution of the BH spin 
direction and magnitude is coupled, with the evolution 
of the spin direction determining the growth or reduction 
of a. 

The paper is organized as follows: in Section 2 we de- 
scribe the gas dynamics and introduce the equations of 

evolution for the mass and spin of the accreting BHs; in 
Section 3 we discuss our results; in Section 4 we summa- 
rize the results and in Section 5 we discuss the potential 
implications of our study on the growth of BHs and on 
jet formation. 

2. METHODOLOGY 

We follow individual BH histories by tracing the evolu- 
tion of the black hole mass Mbh and angular momentum 
vector Jbh = {aGM'^^/c)3-Bn, where we denote with 
< a < 1 the dimensionless spin parameter and with 
Jbh the spin orientation. A single history is a sequence 
of multiple accretion episodes during which Mbh, a and 
Jbh vary upon time t according to the recipes described 
below. 

2.1. The accretion disk properties 

In every single accretion episode the BH is assumed 
to be surrounded by a stationary geometrically thin, op- 
tically thick a— disk (Shakura & Sunyaev 1973). The 
accretion rate M is expressed in terms of the Edding- 
ton factor /Edd, and the accretion efficiency r] (which is 
a fimction of a, Bardeen 1970; Bardeen et al. 1972) ac- 
cording to the relation 



/Edd .^BH 
?7 TS 



= 2x10-^Mbh,6 Moyr-i, 



(1) 

where ?7o.i is the efficiency in units of 0.1, Mbh, 6 the BH 
mass in units of 10^ M©, rs = c7Tc/(47rGmp) = 4.5 x 10^ 
yr, c7t is the Thomson cross-section, and rrip is the proton 
mass. 



We set the total mass of the disk mdisk at each accre- 
tion episode equal to the minimum between (i)mcioud, 
i.e. the mass of a gas cloud that we assume to be available 
in a single feeding episode, and {ii) the self-gravitating 
mass rrisg, the largest for an a— disk to be stable against 
fragmentation by its own self gravity (e.g. Kolykhalov & 
Sunyaev 1980). The mass mdoud is taken to be a con- 
stant, i.e. independent of the BH mass, and carries two 
possible values IO^Mq and 1O'''M0. Within a given ac- 
cretion history, i.e. over the whole BH's life, the mass of 
the cloud is kept fixed at either lO^M© and IO^Mq. The 
self-gravitating mass mgg, for the a— disk, is computed 
from the distance i?disk.sg at which the Toomre parame- 
ter Q ^ CsfJ/TrGE <^ 1; so that at radii R < i?disk.sg the 
disk is stable (where Cg is the central sound speed, 51 the 
Keplerian angular velocity of fluid elements in the disk, 
and S the surface mass density). This yields 



isk,sg 



Rg 



52/45 



/sdd 
VO.l 



-22/45 



(2) 



where ao.i = a/0.1 is the radial shear viscosity param- 
eter in units of 0.1, and Rq the BH gravitational ra- 
dius i?G = "20 Mbh /c'^- Using the solution for the sur- 
face density in the external region of the disk, S(i?) = 

So(i?/i?G)"^/^ with 



(3) 



(4) 



/ f \ 7/10 
-7 in7 -4/5 71^19/20 / /Edd \ _2 

So = 7 X lO'ao./ Mgjle gem ^, 

V ^.i / 

the self-gravitating disk mass reads 

m,g « (87r/5)i?^So(i?disk,sg/i?G)'/', 
and accordingly 

m,,^2x lO^a-'r (^) M^t(f Mo- (5) 
V Vo.i J 

When the disk mass equals mdoud) the outer radius 
-Rdisk.ci is computed as the integral over the surface den- 
sity yielding. 



.Rdisk,cl 

Rg 



4x10* 



^cloud 

lO^M© 



4/5 / f \- 

16/25 ,,-44/25 / /Edd \ 

-0.1 [-^) 
.(6) 

Very massive BHs can stabilize huge accretion disks 
against their own self-gravity. As an example, rrigg ^ 
lO^M© for Mbh ~ 10^ M©. Such a large reservoir of 
mass may not be available in a single accretion episode. 
In fact observations indicate that the mass fraction of 
cold gas relative to stars decreases with increasing galaxy 
mass (e.g. di Scrego Alighieri et al. 2007, CatincUa ct 
al. 2010). To avoid unrealistically massive accretion 
episodes, we adopt the following criterion: if mdoud > 
mgg, the disk around the accreting BH is truncated at 
-Rdisk,sg and has mass mgg; by contrast if mdoud < mgg, 
the disk is truncated at -Rdisk.ci and has mass mdoud- 
Thus for a given value of mdoud, representing the mass 
reservoir in a single accretion episode, there exists a char- 

mdoud ■ 



acteristic BH mass Mgg* for which mgg 



M^S!6 = 0.4ay'' 



^0.1 



.fadd 

Vo.i 



-2/17 



m-doud 



( 104 Mq 



45/34 



(7) 
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Fig. 1. — Distribution of the directions of Jdisk for different degrees of anisotropy. Top-left panel refers to the case of complete isotropy 
("chaotic" accretion, F = 0.5). Top-right, bottom-left, and bottom-right panels refer to accretion with levels of anisotropy F = 0.25, 0.125 
and 0, respectively. 



The mass during a single accretion episode is consumed 
over a timescale Tacc that in the two regimes (when the 
BH mass is below or above Mg^* respectively) reads: 



-11 

BH,6 



/45 / /] 



/Edd 
110.1 



-14/45 



acc,cl 



^cloud 

M 



jEdd 
VO.l 



(8) 

1 

yr- 

(9) 



2.2. The black hole and the disk angular momenta 

The ratio of the angular momentum carried by the 
disk and the angular momentum of the black hole is the 
main driver of how much an accretion episode, and a se- 
ries of accretion episodes can modify the direction and 
magnitude of a BH spin. The a— disk carries an an- 
gular momentum Jdisk obtained integrating over all an- 
nuli the disk angular momentum surface density L(i?) = 
L{R)l{R), where l(i?) is the local orientation of L in the 
disk, and L{R) = j:{R)n{R)R'^ . Following Perego et al. 
(2009; Perego09 hereafter). 



Jdisk(i?) oc M^/GM^R''/\ 



(10) 
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The ratio between the angular momentum of the disk 
and BH reads; 



Jdisk(^disk,sg) „ „ 13/45 / /Edd 
Jbh \V0.1 



-7/45 



,,r-37/45 -1 



(11) 

Since Jdisk c< i?^/^, disks that are not truncated by seh- 
gravity carry an angular momentum smaller by a factor 
(^disk,ci/^disk,sg)^^'* than disks with mass TOsg, and this 
will play a role in modeling the BH history. 

The orientation of Jdisk, denoted as Jdisk, is deter- 
mined mainly by the angular momentum in the outer- 
most region of the disk. Jdisk is uncorrelated with the BH 
spin direction Jbh at the onset of any accretion episode. 
Hereon, disk orientation is described in terms of the polar 
angle ^disk between Jdisk and a fixed reference direction, 
i.e. the unit vector in our spherical coordinate refer- 
ence system: 6'disk = cos" ^ (Jdisk ■ e^). 

The orientation of the accretion disk is selected 
through Montecarlo sampling. We follow four differ- 
ent prescriptions that correspond to different degrees of 
anisotropy in the BH fuelling process. In the first, Jdisk 
is distributed at random. This corresponds to a uni- 
form distribution for the azimuthal angle 0disk and to 
a distribution proportional to sin 9 for the polar angle 
^disk- In this particular case, the fraction F of accre- 
tion events with ^disk > 90° is 0.5. The other three 
cases are constructed following initially the same pro- 
cedure, but inverting at random the sign of ^disk to 
build three distributions with F = 0.25, F = 0.125, 
and F = 0. F = 0.25 corresponds to a distribution 
with three times more events in the "northern hemi- 
sphere" (0disk < 90°) than in the "southern hemisphere" 
(^disk > 90°), F = 0.125 corresponds to a distribution 
with a north-to-south events ratio of seven, and the dis- 
tribution with F — has no events in the southern hemi- 
sphere. Figure [1] shows the distributions of Jdisk for the 
four values of F considered. Note that the F = 0.5 case 
is exactly what is assumed in the standard "chaotic ac- 
cretion" scenario. F = 0, on the other hand, does not 
correspond to "coherent accretion" , but it mimics accre- 
tion through disk clouds that are distributed isotropi- 
cally but that share a common sense of rotation. F — 
corresponds to a 3D-dispersion to rotation velocity ratio 
f/wrot 1 for the gas fuelling the BH. 

The direction of the BH spin vector Jbh , initially se- 
lected at random, is followed by tracking the BH spin and 
accretion history and by recording the values of a, Jbh 
and Mbh at the end of each accretion episode. Hereon 
we denote with 0bh the polar angle of Jbh with e^, i.e. 
^bh = cos" ^ (Jbh • e^). The BH spin vector can point in 
all directions, and Jbh as well referred to the 

reference frame of the galaxy defined by e^. 

2.3. a, Jbh and A/bh on a history path 

At the onset of any accretion episode the BH spin 
vector Jbh is generally misaligned with respect to the 
direction of the accretion disk Jdisk- In high viscosity 
a— disks, this configuration is unstable and evolving into 
a lower energy state. There is an early phase in which 
gravito-magnetic torques exerted by the spinning BH on 



disk's fluid elements cause the disk to warp on a timescale 
T'warp (Bardeen & Petterson 1975); the second phase is 
the alignment phase, i.e. a change in the orientation of 
the BH spin over a longer time. 

In the early phase, the spinning BH induces Lense- 
Thirring precession of the orbital plane of disk's fluid el- 
ements: precession of the plane occurs at a frequency 
cjp.ec = (2G/c2)JBH/i?^ = {2G/c)aMl^/R^, so that 
fluid elements closer to the BH precess faster. Close to 
the BH, their orbital plane tends to align (or anti-align if 

counter- rotating) parallel to the BH spin direction Jbh, 
on times shorter at shorter radii, so that the perturbation 
diffuses radially outwards, on a timescale r^arp ^ /'^2, 
where 1^2 is the vertical shear viscosity. The disk is max- 
imally warped at i?warp, corresponding to the distance 
where the warp-diffusion timescale becomes comparable 
or shorter than the precession time w~4c- 



4GJ 



BH 



Rg 



V2C^Rg 



500a 



24/35^4/7^^4/35 



0.1 



'BH,6 



/Edd 
?70.1 



,4/7 



(12) 

where V2 ^ a2CsH with a2 ~ jy^ji^a) (Lodato & Pringle 
2007; PeregoOQ for details) i?warp indicates the region in 
the inner disk where the fluid elements align or antialign 
with Jbh- kX. R 

warp the warp timescale is 



' warp 



„(. 72/35 .-12/7.^47/35 / /Edd 



5/7 

a ' yr. 



. (13) 

Since the warp timescale is shorter than the viscous 
timescale at all annuli, the deformation diffuses more 
rapidly outwards than the inward radial drift motion so 
that the deformed disk attains an equilibrium profile, and 
the shape of the perturbation is stationary (Martin et al. 
2007). 

The warped disk (in the small deformation approxi- 
mation) is described by an equilibrium surface density 
E(i?) and a velocity (l{R) close to that of an unper- 
turbed Keplerian disk, plus a deformation in the lo- 
cal angular momentum vector (per unit surface area) 

L = L{R){lx,ly,lz) (with 1 a unit vector in the direc- 
tion of L; we defer to Perego09 for details). Angular 
momentum conservation then imposes the spinning BH 
to precess and align relative to the total angular mo- 
mentum Jtot = Jdisk + Jbh- During individual accretion 
episodes at a rate M, the BH mass increases according 
to: 



= (1 - r^)Mc^ = {1-V)( (14) 



dt 

with an e— folding timescale 



5x10'- yr. (15) 

(l-7?)./Edd (1 - fyo.lj/Edd 

The BH angular momentum vector evolves according to 

L X Jbh 



—7—=^ Aiscol(i?isco)+^- / 55 — dR 

at C C Jdisk ^ 

(16) 

where the first term in equation (15) describes the change 
in the spin modulus a due to the transfer of angular 
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momentum per unit mass at the innermost stable circular 
orbit i?isco (where A(i?isco) is a known dimensionless 
function of a). Gravito-magnetic coupling in the inner 
region of the precessing disk ensures that the direction 
of l(i?isco) is parallel o anti-parallel to Jbh- 

The second term of equation describes the change 
in the BH spin orientation which tends to reduce the de- 
gree of misalignment between the disk and the BH spin 
(Perego09; Scheuer & Feiler 1996). An intuitive expla- 
nation for the alignment is as follows. The warping of 
the disk results always in a either aligned or antialigned 
inner part of the accretion disk, due to the BH metric. 
On the other hand, since the frequency of the Lense- 
Thirring precession decreases with radius, the outer part 
of the accretion disk can effectively be considered un- 
perturbed (this corresponds to the Newtonian limit for 
the accretion disk), and there exists a region (near the 
so called warp radius) where the local disk angular mo- 
mentum unit vector is strongly misaligned both with the 
spin and disk's angular momentum at far distances (in 
the unperturbed disk). In other words, gas moving from 
the unperturbed outer region into the innermost region 
modifies its angular momentum direction. As a conse- 
quence, the BH has to (at least partially) align with the 
direction of the angular momentum of the outer accretion 
disk to ensure conservation of the total angular momen- 
tum, and to compensate the "loss/change" of angular 
momentum in the warp region. The alignment effect on 
the BH spin due to the cumulative torques exerted by a 
mass distribution does not necessarily require the pres- 
ence of viscosity. The same alignment can be induced 
(under some conditions) by the spin orbit interaction be- 
tween the BH and a rotating stellar cusp (see Merritt & 
Vasiliev 2012). Viscosity guarantees alignment of the in- 
ner disk region and thus guarantees the occurrence of a 
stationary warp region in the disk causing the torquing of 
the BH and in turn spin alignment. From equation (|16p . 
we can infer an evolution equation for the spin modulus 



da 
'dt 



A{R 



ISCOJ 



2a 



1 



1 



Edd 



(17) 



Equation ()17p shows that changes in the spin modulus 
a typically occur on the timescale Tgpin = a/ a compa- 
rable to tmbh (^1- HSI- As shown in Bardeen (1970), a 
non-spinning BH (a = 0) is spun up to an extreme Kerr 
BH (a = 1) after having accreted a mass V6AfBH- The 
timescale of the alignment process (jl6[) is much shorter. 
The peak of the torque perpendicular to the BH spin 
(that responsible for the evolution of the spin direction) 
is maximal at the warp radius (where the direction of the 
angular momentum of the gas in the disk changes). The 
warp radius is considerably larger than the last stable or- 
bit, so the angular momentum per unit of mass of the gas 
at the warp radius is much larger than that at the ISCO. 
As a consequence, the torque exerted onto the BH re- 
sponsible for the evolution of its direction is considerably 
larger than that responsible for the spin magnitude evo- 
lution. Consequently, this results in a shorter timescale 
for the spin alignment with respect to that associated to 
the spin magnitude evolution. 

Equations (fT4|) and (|16|) are integrated considering ini- 
tial BH masses of A/bh,o = 10^ M0, and arbitrary ini- 
tial spin moduli (oq) and orientations; /Edd that en- 



ters all timescales can be considered as a scaling pa- 
rameter so that smaller values of /sdd imply longer 
times. In the calculation we fix /Edd = 0.1. Accord- 
ing to equation (|16|) , if the disk angular momentum con- 
tributes most to Jtot, the BH spin-disk orientation an- 
gle CsH.disk = cos~^(Jbh ■ Jdisk) reduces to zero (i.e. 
full alignment even starting with CfiH.disk = 180°) on 
a timescale 



Tal 10 ao,i A'^BH,6 



35 



/Edd 
'70.1 



a^/V (18) 



(Martin et al. 2007; Lodato & Pringle 2006; Perego09). 
The timescales r^arp and Tai depend on a and A/bh ex- 
plicitly, and the following inequalities hold 



spin ■ 



(19) 

In our study, we integrate the BH spin evolution equation 
over timesteps 5t that are r^arp < 5t < Tai tracing con- 
temporarily the change in a, Jbh and A/bh- We notice 
that as Tspin > Tai, the spin orientation changes rapidly 
with time, so that we can exclude the possibility that ret- 
rograde accretion reduces a to before re-increasing it to 
1. The rapid change in Jbh compared to the change in 
a, favors conditions of prograde accretion episodes, even 
starting from retrograde conditions. 

The progressive increase in A/bh during accretion 
episodes with constant Wdisk = "icioud (when A/bh > 
A/gy ) implies the formation of an accretion disk with 
'7^disk/A/BH always decreasing. Thus the disks that form 
during single accretion episodes are progressively smaller 
and carry less angular momentum relative to the BH. Un- 
der these circumstances and depending on a, the dimen- 
sionless ratio Jdisk/ Jbh ^ 1, and the gravito-magnetic 
BH-disk coupling has little influence on the BH spin di- 
rection. At even larger masses (already in the regime 
where Jdisk/Jen < 1) the warp radius /?warp may rise 
above > /?disk depending on a, and equation be- 
comes invalid (Martin et al. 2007). This transition oc- 
curs when A/bh, 6 > -^bh'^6' where 



71 rwarp 1 n I 'Ttcioud 



35/82 



-1/41-25/82 



^0.1 



/Edd 



.(20) 

Under these circumstances, we assume that in any single 
event Jbh aligns instantaneously to the direction of the 
total angular momentum Jtot (King et al. 2005). The 
re-orientation of Jbh is rather small, since in this regime 
the BH spin is quite close to the total angular momen- 
tum already at the beginning of every individual accre- 
tion episode (Jtot is in fact dominated by Jbh)- Jdisk 
instead undergoes a fast and significant re-orientation, 
and is either aligned or antialigned with respect to Jbh 
depending on the initial angle of relative misalignment 
(Cnn.disk) and the Jdisk/ Jbh ratio. The two angular mo- 
menta are aligned if cos(CBH,disk) > -Jdisk/2JBH (or 
counter- aligned when the opposite inequality holds; King 
et al. 2005; Lodato & Pringle 2006). The BH mass 
and spin parameter are then evolved according to equa- 
tion (HI]) and P7)) . 

To summarize, our procedure leads to the following 
sequence: 

• The disk mass is the minimum between the cloud mass 
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and the self-gravitating mass, mdisk = niin(mcioud, ?7isg)- 
The transition between the two phases occurs when the 
BH mass reaches the value M^^q given by equation ([7]); 
then the disk contains the whole cloud mass. 
• Up until Mbhs < M^^Q, we follow the spin-disk align- 
ment as a function of time by solving for equation (IT6)) . 
Here, M^^^ is given by equation ([20| . Above this mass 
the disk is assumed to be aligned or anti-aligned with the 
angular momentum of the BH, depending on the initial 
angle of relative misalignment. 
Therefore, three regimes exist, Mbhs < M^^^; 

M^K,6 < ^BH,6 < M^^J, and Mbh,6 > A^eS'e- While 
the switch between the first and second regime depends 
only on the properties of the gas fuelling the BH, the 
switch between the second and third regime depends also 
on the BH spin, and this threshold mass, therefore, de- 
pends on the whole growth history of the BH, i.e., how it 
gained its spin. We anticipate that, regardless of the 
regime, the important parameter that determines BH 
spin evolution is the ratio Jdisk/JBH, as previously dis- 
cussed. 

3. RESULTS 

Figure [2] shows the evolution of the ratio between Jdisk 
and Jbh as a function of the BH mass, for accretion 
histories that differ in the value of rricioud and of the 
anisotropy parameter F. 

As discussed in the previous section we can distin- 
guish an early phase (I) during which the accretion disks 
carry a large angular momentum, and a second phase 
(H) where the opposite holds. During phase I, the disks 
are initially truncated by their own self-gravity and their 
mass is determined by the BH mass. As the BH grows 
in mass the disks carry the cloud mass mdoud- This 
transition occurs at and is visible as a knee in the 

Jdisk/ Jbh versus BH mass diagram. During phase H, 
disks are tiny (-Rwarp > -Rdisk) and there is the switch 
between the two prescriptions for the BH spin evolution. 
This occurs around a few 10^ Mq and is highlighted as a 
red shaded area in each panel. The upper left and right 
panels refer to the isotropic case {F — 0.5) with a max- 
imum mass per accretion event Wcioud = 10^ Mq and 
10^ Mq, respectively. The minimum value of Jdisk/ Jbh 
for the two cases is between 10~^ — 10~^ and is attained 
when Mbh ^ 10^ M©. Similar results are found for differ- 
ent values of the anisotropy parameter, as shown in the 
bottom panels for mdoud = 10^ Mq, assuming F — 0.25 
and respectively. Notice that in these cases Jdisk/ Jbh 
drops below 10~^ for very large masses and this reflects 
the large spin carried by these BHs. For anisotropic cases 
{F = 0.25; 0) the massive BHs tend to align their spins 
with the average angular momentum of the gas reservoir, 
as will be discuss in section [SHJ This results in an higher 
fraction of prograde accretion events, and, as a conse- 
quence, in a higher Jbh- A detailed description of the 
evolution in direction and magnitude of a is presented in 
the next two sections. 

3.1. Evolution of the BH spin direction 

The BH spin orientation in a fixed reference frame is 
defined by two angles, the polar angle ^bh, and the az- 
imuthal angle (pBH- Figure [3] shows, as an example, the 



evolution of 0bh as a function of the BH mass for differ- 
ent accretion histories with F — 0.5 (upper panels) 0.25 
and (lower left and right panels respectively), and cloud 
masses mdoud = IO^'Mq (upper left panel), and IO^Mq 
(other panels). In case of perfect isotropy in the BH fu- 
eling {F = 0.5; upper panels in Figure [2]), the lack of 
any preferential direction results in an unbiased random 
walk of the BH spin direction. This can be quantita- 
tively viewed by plotting the distributions of 9bb and 
0BH averaged over 500 accretion histories extracted at 
random. In the upper panels of Figure |4] we show the 
distribution of the two angles for the isotropic case for 
BHs in the mass interval 10^ Mq — 10^ Mq (green, long 
dot-dashed line), 10^ Mq - lO"^ Mq (purple, short dot- 
dashed line), 10^ Mq - 10* Mq (red, long dashed line), 
10* Mq - lO'^ Mq (blue, short dashed line) and > IO^Mq 
(black, solid line). It is clear that both distributions are 
consistent with isotropy. 

The lower-left panel of Figure [3] refers to the case 
when F = 0.25, corresponding to 3/4 of the accretion 
episodes having the disk angular momentum confined 
in the northern hemisphere. The behavior of 0bh dif- 
fers in the two accretion phases (described in the previ- 
ous section). When Jdisk /Jbh > 1 (for masses Mbh ^ 
few 10^ Mq) the spin orientation aligns significantly with 
the angular momentum of each accretion episode, thus 
changing erratically over the full range (0, 180°, despite 
the fact that the mean is less than 90°). By contrast, dur- 
ing phase II (Jdisk/JBH < 1) the spin tends to align with 
the average angular momentum of the accreting material. 
In this phase the BH spin slightly tilts toward the direc- 
tion of the angular momentum of each accreting cloud. 
The sequence of accretion episodes results in a random 
walk biased toward the direction of the average angu- 
lar momentum of the inflow. This trend is illustrated in 
the two middle panels of Figure ID While 0bh retains a 
flat distribution, the distribution of ^bh covers the full 
[0, 180°] range for small BH masses, with a 3 times higher 
probability of 9bh < 90°. For higher masses (i.e. dur- 
ing phase II) , the distribution of 9bh shifts toward lower 
angles, being confined below ^ 60° for Mbh ^ IO^Mq. 
Similarly, for the most anisotropic case considered here 
{F = O), the BH spin is confined within a solid angle 
that decreases with increasing BH mass (see the lower 
right panel of Figure [3] and the lower two panels of Fig- 
ure H]). For the very massive BHs (Mbh > 10^ Mq), Jbh 
is nearly constant, demonstrating the stability of the spin 
direction even during substantial increase in mass (from 
10^ to 5 X IO^Mq). 

In summary, we show that for small BH masses 
(Mbh ^ 10^ Mq) the BH spin always aligns to the disk 
angular momentum in every single accretion episode, and 
can be substantially misaligned relatively to the aver- 
age of the angular momenta of the disks. For Mbh ^ 
10'' Mq, a single accretion episode does not modify signif- 
icantly the BH spin direction. In this regime, the BH spin 
aligns with the average direction of the angular momen- 
tum of the accreting material, with a degree of alignment 
that increases with increasing anisotropy in the fueling 
process. 

3.2. Evolution of the BH spin magnitude a 
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Fig. 2. — Ratio between the total angular momentum of the accretion disk Jdisk Eind the BH spin Jbhi a function of the BH mass 
Mbh I along an accretion history. The upper left panel refers to the isotropic case (F = 0.5) with a maximum mass per accretion event 



^cloud 

respectively. Three regimes exist, Mbh, 6 < -'^^bh* 6 



IO^Mq and F = 0.5, 0.25, and 0, 



The results in the upper right, lower left and lower right panels assume mdoud 

■jij'g < Mbh, 6 < ^bh'^6' ^'^'^ AfBH,6 > -'^bh'^6' "^^^ knee at low BH masses shows 
the transition between the first and the second regime. The second transition is highlighted by the shaded vertical area (see text for details). 



In this section we focus on the evolution of the spin 
magnitude a considering for the first time its coupling 
with the orientation Jbh and with the dynamical prop- 
erties of the accretion. 

Figure [5] shows the evolution of the BH spin a as a 
function of its mass, for the isotropic case and two dif- 
ferent values of mcioud- The memory of the initial spin 
is erased after the BH accretes a few times its initial 
mass (AfBH,o — 10^ Mq), and a rises attaining values 
~ 0.9 during phase I. The large spin is a consequence of 
the rapid alignment induced by the Bardeen-Petterson 



effect that occurs on a timescale Tai < Tacc < ^spin 
turning initially retrograde accretion into prograde accre- 
tion before disk consumption. Later, i.e. for larger BH 
masses, the spin drops to lower values, down to a ~ 0, for 
^BH ^ 10^ Mq. The beginning of the spin decrease cor- 
responds to the transition between phase I and H that 
occurs at different BH masses, depending on the value 
of mdoud, as shown in the figure. The parameter that 
controls this transition is Jdisk/JsH: smaller clouds carry 
lower angular momenta and the transition between phase 
I and H occurs at smaller BH masses. The spin mod- 
ulus decreases because over a single accretion episode 
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Fig. 3. — Evolution of the polar angle Bbh describing the BH spin orientation, relative to a fixed coordinate system, as a function of the 
BH mass. The upper left panel refers to the isotropic case {F = 0.5) with a maximum mass per accretion event rricloud = IO^Mq. The 
results in the upper right, lower left and lower right panels assume rrtdoud = IO^Mq and F = 0.5, 0.25, and 0, respectively. 



Jbh does not change significantly. Retrograde accretion 
episodes remain retrograde over the disk consumption 
timescale, so that, for isotropic fueling, the probability 
of having a prograde accretion events is exactly the same 
of having a retrograde one. Because the location of the 
IS CO is farther away, and therefore the accreted mate- 
rial carries a larger specific angular momentum, retro- 
grade accretion transfers more angular momentum per 
unit mass than prograde accretion, isotropic fueling re- 
sults in net spin-down and thus low spins (e.g. King & 
Pringle 2006). 

Figure [6] shows the spin evolution for different value of 
the anisotropy parameter F — 0.5; 0.25; 0.125; (with 



Ticioud — 10^ Mq). The early phase (I) is similar in all 
the cases. Rapid BH-disk alignment results in a « 0.9. 
At the transition between phase I and II, the spin always 
decays, because at this transition the BH-disk alignment 
is only partial, and retrograde accretion is possible for a 
significant fraction of the time that elapses during each 
accretion event. By contrast, the value of the BH spin 
at large masses {Mbh ^ 10^ M©) is strongly dependent 
on the degree of anisotropy. For F = 0.25, large BHs 
carry spins with a ~ 0.45, and this asymptotic value of a 
increases with decreasing F. For F = 0, i.e. the highest 
degree of anisotropy described here, a — > 1. This is a con- 
sequence of the stability of the spin direction (discussed 
in section [3?T|) . For highly anisotropic cases, the massive 
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Fig. 4. — Distributions of the polar angle 5bh (left column) and azimuth angle (/)bh (right column), for BHs of different masses, assuming 
"^cloud = 10^ M0. Top, middle, and bottom panels refer to the F = 0.5, 0.25, and cases, respectively. BHs in the mass interval 
10^ Mq - 10*^ M0 are denoted with green, long dot-dashed line, 10'' Mq - lO'' Mq with purple, short dot-dashed line, lO'' Mq - 10* Mq 
with red, long-dashed line, 10* Mq — 10^ Mq with blue, short-dashed line, and > 10^ Mq with black, solid line. 



BHs tend to align their spins with the average angular 
momentum of the gas reservoir, resulting in an higher 
fraction of prograde accretion events. Thus, we can cor- 
relate the spin of the large BHs with the dynamics of the 
fueling mechanism. 

Figure [7] shows the distribution of the BH spin in dif- 
ferent mass intervals, assuming F = Q and mdoud — 
10^ M0. Here we emphasize that the highest spins are 
attained at the largest BH masses. The distribution 
corresponding to 10^ M© < Mbh < 10^-^ Mq peaks at 
a ~ 0.997. This may have important implications on 
the efficiency at which these BHs can accrete, on the 



10 

statistics of very massive BHs, and on the possibility of 
launching jets (e.g. Tchekhovskoy & McKinney 2012). 
A large degree of anisotropy {F = 0) is the only way to 
get maximally rotating BHs at the high mass end. 

4. COMPARISON WITH OBSERVATIONS 

It is not straightforward to compare the predictions 
of our simple investigation with constraints on the BH 
masses and spins inferred from observations. Theoreti- 
cally, a clear prediction of how the expected spin distribu- 
tion evolves with the BH mass can be made once the dis- 
tribution/dynamics of the fueling process is known (i.e., 
when F and mdoud are known, in our context), and fuel- 
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Fig. 5. — Evolution of the BH spin magnitude a as a function of 
-'^BHi for the isotropic case (_F = 0.5). The upper (lower) panel 
refers to accretion episodes with rricloud = IO^Mq (10'*Mq). The 
black line refers to the mean over 500 realizations. Red and or- 
ange shaded areas enclose intervals at 1 — tr and 2 — a deviations, 
respectively. 

ing can in principle differ depending on the galaxy types 
and masses. Observationally, the constraints on BH spins 
are often debated. However, we can make some very 
general comments on how our model fits within the ob- 
served trends. As discussed in the Introduction, the spin 
of a number of massive BHs has been measured through 
spectral fitting of the broad Ka iron lines at 6.4 keV. In 
almost all the cases, the measured spins are larger than 
0.5, for BHs of masses 10^ Mq < Mbh ^ 10* Mq. The 
spin parameters of the BHs hosted in Swift J2127, NGC 
3783, and MGC-6-30-15 are still debated, and could be 
lower than 0.5, with a minimum of a < 0.32 for NGC 
3783 (Patrick et al. 2011b, but see also Brenneman et 
al. 2011 for a higher estimate of a). These measurements 
are in agreement with our simple model, i. e. there is al- 
ways a history that can recover the observed spin for the 
mass observed. In particular these data suggest a non- 
negligible anisotropy in the fuelling of these BHs. Note 
however that this could be related to an observational 
bias toward measurements of highly spinning objects, as 
discussed in Brenneman et al. (2011). 

A possible evolution of the BH spin with mass is found 
in the countinuity-equation based study by Shankar, 
Weinberg & Miralda-Escude (2011). In this study, a ra- 
diative efficiency (and as a consequence, a spin param- 
eter) increasing with the BH mass is needed to avoid 
overproduction of BHs at the high mass end of the local 
observed mass function. A similar study by Li, Wang & 
Ho (2012) finds i] increasing for increasing A/bh at high 
redshift (z > 1). We notice, however, that Li et al. do 
not recover this trend at low redshift. They found the 
radiative efficiency at z < 0.8 to be almost independent 
on the BH mass. The accretion efficiency could be de- 
rived fitting observational data with accretion disk mod- 
els for individual sources (Davis & Laor 2011; Raimundo 



et al. 2012). However, the values determined are largely 
sensitive to the unknown parameters and uncertainties, 
which makes it very hard to draw any conclusions on the 
spin dependence with black hole mass, as discussed in 
Raimundo et al. (2012). 

If the trend of r\ suggested by Shankar et al. (and by Li 
et al. for high redshifts) is confirmed, at least some mas- 
sive BHs (Mbh ^ lO^Mo) must be very close to maximal 
rotation (a > 0.98), because of the very steep evolution 
of r\ close to a « 1. The dependence of on a is quite 
sensitive as r\ varies from 0.151 for a = 0.90, to 0.43 for 
a = 1. BHs can be spun up to such high spins (a > 0.98) 
through BH-BH binary coalescences only when specific 
conditions on the spins of the parent BHs and the orbital 
configuration are met. A coalescing binary can result in 
a remnant BH with awl only if i the parent BHs were 
close to maximal rotation prior to merge, and had their 
spins aligned with the binary orbital angular momen- 
tum (Marronetti et al. 2008; Kesden et al. 2010al], ii 
through a long sequence of extremely unequal mergers 
all confined in the same well defined orbital plane (cf the 
"equatorial case" in Berti & Volonteri 2008). For coalesc- 
ing BHs not rapidly rotating or with spins not aligned 
with their orbital angular momentum, the remnant is ex- 
pected to have a smaller spin, with an average of a ^ 0.7 
(Hughes & Blandford 2003; Berti & Volonteri 2008). To 
summarize, BH coalescences can preserve high spins, but 
cannot be the original cause of extreme spins (Berti & 
Volonteri 2008). The reason for close-to- maximally ro- 
tating BHs must lie in their accretion history. The very 
large spins in the most massive BHs, if confirmed, will 
be an indication of a net degree of anisotropy in the fu- 
eling of this objects (F w 0). We note that our analysis 
predicts for the first time that the most rapidly rotating 
BHs (a > 0.98) are also the most massive ones, and in 
addition, that the spin orientation is stable over many 
accretion cycles. 

An independent constraint comes from jet observa- 
tions. The bulk of luminous radio-loud AGN is associ- 
ated with very massive BHs (e.g. McLure & Jarvis 2004; 
Metcalf & Magliocchetti 2006; Shankar et al. 2008a,b; 
Shankar et al. 2010) typically hosted in very massive 
ellipticals (e.g. Capetti & Balmaverde 2006). If the 
radio jet is powered by the extraction of rotational en- 
ergy of the BHs, powerful radio sources would be re- 
lated to rapidly rotating BHs (e.g. Blandford & Zna- 
jek 1977). Recent MHD simulations suggest that the 
power is a steep increasing function of the BH spin for 
a ~ 1 (e.g. Tchekhovskoy & McKinney 2012, and refer- 
ences therein). In this scenario, our model predicts the 
most massive BHs to produce the most powerful jets, 
and jets are stable in their orientation if the BHs are 
not completely isotropically fueled (but not necessarily 
experience coherent accretion) . The high mass BHs that 
have been fed by almost isotropically distributed gas have 
lower spins, and lower radio power. The difference in the 
fuelling process could explain the observed dichotomy of 
the radio-loudncsses of AGN (Sikora, Stawarz & Lasota 
2007). The prediction of the most powerful jets being re- 
lated to the most massive BHs is also in agreement with 

^ This configuration is predicted for gas rich galaxy mergers 
(Bogdanovic, Reynolds k. Miller 2007; Dotti et al. 2010; Kesden 
et al. 2010b; Berti, Kesden & Sperhake 2012). 
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Fig. 6. — Evolution of the BH spin magnitude a as a function of Mbhi assuming rridoud = IO^Mq. The upper left, upper right, lower 
left, and lower right panels refer to _F = 0.5, 0.25, 0.125, and 0, respectively. The colour code is as in figureO 



the observed trend of the increasing fraction of radio- loud 
AGN for increasing BH masses (McLure et al. 1999; 
Jiang et al. 2007; Caccianiga et al. 2010; Chiaberge 
& Marconi 2011), and with the evidence that the most 
powerful radio and 7— loud AGN are associated to the 
heaviest supermassive BHs in the universe, with masses 
> lO^M© (Ghisellini et al. 2010a,b). 

We can also compare the distributions of the spin ori- 
entations for BHs of different masses as predicted by 
our model, with the observational constraints available 
to date. It is commonly accepted that the directions of 
relativistic jets in radio-loud AGN are not in a close cor- 
relation with the galaxy morphology (e.g. Schmitt et al. 
2002; Verdoes Kleijn & de Zeeuw 2005). However, dif- 



ferent degrees of alignment are observed, depending on 
the morphologies (and masses) of the host galaxies. Re- 
garding Seyfert galaxies, the relative orientation of the 
jets and the minor axes of the galaxy is consistent with 
being drawn from an isotropic distribution (Clarke, Kin- 
ney & Pringle 1998; Nagar & Wilson 1999; Kinney et 
al. 2000). The degree of alignment seem to increase 
moving toward more massive galaxies and more massive 
BHs. Schmitt et al. (2002) studied the degree of mis- 
alignment between the direction of jets in radio-galaxies 
(early type and SO) and the direction normal to their nu- 
clear disks, as traced by dust lanes. They found that the 
jets are typically misaligned with respect to the normal 
to the nuclear disks by less than 55° ~ 80°. Similar re- 
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Fig. 7. — Distribution of the BH spin magnitude a for BHs in six different mass intervals, for F = and mdoud = IO^Mq. 



suits have been obtained by Verdoes Kleijn & de Zeeuw 
(2005) for low-power radio galaxies. They found slightly 
better alignments when comparing the direction of the 
jet and the normal to ellipsoidal-shaped dust lanes (trac- 
ers of nuclear disks), with a peak in the distribution of 
misalignments at ~ 45°. Such a level of alignment is in 
contrast with an isotropic distribution, where a peak of 
misalignments is expected to coincide with the peak of 
the solid angle, at ~ 90°. Note however that, given the 
uncertainties in the estimates, an isotropic distribution 
cannot be ruled out at more than the 95% confidence 
level (Verdoes Kleijn & de Zeeuw 2005). Verdoes Kleijn 
& de Zeeuw suggest that while dust ellipses trace gas on 
unperturbed Keplerian orbits that does not contribute 
to the fueling of the central AGN, dust lanes can bet- 



ter trace perturbed gas that traces the material falling 
toward the BH. The degree of alignment between the 
normal to dust lanes and the jet direction is sensibly 
higher, with a peak of the distribution at 20° ^ 30°. 
These studies suggest that a relation between jets and 
morphological properties of the galaxies exists only for 
massive galaxies, hosting massive BHs. This result is in 
broad agreement with our model (see Figure [J]) that also 
indicates the stability of the jet direction, stability that 
is necessary to launch a large-scale jet, as seen in radio 
galaxies. We notice that on average the direction of the 
jet is expected to be related to the morphology of the nu- 
clear distribution of gas, and not necessarily to the large 
scale stellar distribution. However, Saripalli & Subrah- 
manyan (2009) found that the largest radio sources (with 
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a projected linear size exceeding 700 kpc) have the jets 
preferentially aligned with the minor axes of their host 
galaxies. In the framework of our model, such a sta- 
ble direction of very large scale jets can be explained if 
the BH is sufficiently massive (> lO^M©), and hosted 
in a galaxy with a significant degree of anisotropy in its 
gas component. This trend is observed also by Browne 
& Battye (2010), as they observe ellipticals with a net 
rotation hosting jets preferentially aligned with the ro- 
tation axes, while non-rotating/triaxial ellipticals do not 
show such an alignment. 

The predicted dependence of the magnitude and di- 
rection of Jbh on F can be further tested studying the 
nuclear gas dynamics around BHs with either an esti- 
mate of the spin parameter or a well constrained jet di- 
rection. Owing to its extreme angular resolution, the 
Atacama Large Millimeter/submillimeter Array is the 
likely to be instrumental in constraining the dynamics 
of dense molecular gas (and consequently F) in the nu- 
clei of nearby galaxies. 



5. SUMMARY AND DISCUSSION 

In this paper we explored the evolution of the spin 
parameter a (or Jbh) of massive BHs considering the 
contemporary evolution of the orientation Jbh, along se- 
quences (histories) of accretion episodes. Histories are 
modeled as a succession of single accretion events where 
an a— disk of given mass mdisk and orientation Jdisk 
forms. Disk orientation is not fixed but it is drawn from 
a distribution that carries a degree of anisotropy. This 
anisotropy is expected to be seeded in the gas clouds that 
surround the massive BH in the galactic nucleus and that 
are accreted. Our findings and their astrophysical con- 
sequences can be summarized as follows: 

• At low BH masses (Mbh ^ 10^ Mg) the BH spin al- 
ways aligns to the disk angular momentum in every single 
accretion episode, and can be substantially misaligned 
relatively to the average of the angular momenta of the 
disks. For Mbh ^ 10'' M©, a single accretion episode 
does not modify significantly the BH spin direction. In 
this regime, the BH spin aligns with the average direc- 
tion of the angular momentum of the accreting material, 
with a degree of alignment that increases with increasing 
anisotropy in the fueling process. 

• Small BHs (Mbh ^ 10^ M0) carry large spins a ^ 0.9. 
The spin direction changes erratically from episode to 
episode so that the vector Jbh exhibits a random walk 
behavior, regardless the properties of the fueling process. 

• Large BHs (A/bh ^ 10^ Mq) can carry either low spins 
or large spins, depending on the fueling conditions. The 
spin is low, a sa 0, if the distribution of clouds in the 
hosts is completely random and isotropic. The spin is 
larger, up to a > 0.99, and the spin orientation is stable 
on the sky if the gas accreting onto the central BH posses 
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some degree of anisotropy, i.e. ij the accreting material 
has, on average, non zero angular momentum. Only if 
the degree of anisotropy is high [F = 0) the most massive 
black holes can be maximally rotating. 

• The most anisotropic case [F = 0) studied here does 
not correspond to what is usually referred to as coherent 
accretion: for F = 0, two successive accretion events can 
be misaligned by up to 180 degrees. F = mimics ac- 
cretion through disk clouds that are distributed isotrop- 
ically but that share a common sense of rotation. F — Q 
corresponds to a 3D-dispersion to rotation velocity ratio 
c/wrot ^ 1 for the gas fuelling the BH. 

• There may exist an interval of BH masses, in between 
the two populations described above, in which the spins 
can be far from either zero or unity and changes orienta- 
tion at random. 

• Although our analysis predicts that most low mass BHs 
have substantial spins (a ~ 0.9; see the two bottom pan- 
els of Figure!?]), the most rapidly rotating BHs (a > 0.98) 
are also the most massive ones (see upper panel of Fig- 
ure [7]). In addition, the spin orientation of the most 
massive BHs remains stable over many accretion cycles. 

• The very large spins in the most massive BHs, if con- 
firmed, will be an indication of a net degree of anisotropy 
in the fueling of this objects (F « 0). In other words, our 
model predicts that a clear correlation exists between the 
kinematics of clouds feeding a BH and its spin, but only 
in the case of the most massive BHs (A/bh <L 10^ Mq). 

As a final comment, we notice that light BHs (Mbh ^ 
10^ Mq) carry large spins undergoing erratic changes in 
their orientation. This can have strong implications on 
the efficiency of feedback exerted by active BHs onto 
their host galaxy, that can potentially set the BH-host 
galaxy scale relations (see Nayakshin, Power & King 
2012 and references therein). The higher "inertia" of 
the most massive BHs, if embedded in anisotropically 
moving gas, could reduce the feedback efficiency, since 
only a small solid angle would be affected. By contrast, 
if any anisotropic feedback (such as jets or biconical out- 
flows, etc.) are launched during single episodes around 
less massive BHs, the spin random walk would result in 
a spread of the injected energy over 47r during their life- 
time. 

Future studies will address the consequences of the 
model in the context of the cosmological evolution of 
BHs. 
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